Robotics  Research 
l^hnical  Report 


■■■■^'' 


Polynomial  Remainder  Sequences 
and  Theory  of  Subresultants 

by 

Chung-jen  Ho 
Chee  Keng  Yap 

Technical  Report  No.  319 

Robotics  Report  No.  119 

September,  1987 


\ 


New  York  University 
nt  Institute  of  Mathematical  Sciences 

Computer  Science  Division 

25 1  Mercer  Street   New  York,  N.Y  1 00 1 2 


Polynomial  Remainder  Sequences 
and  Theory  of  Subresultants 

by 

Chung-jen  Ho 
Chee  Keng  Yap 

Technical  Report  No.  319 

Robotics  Report  No.  119 

September,  1987 


New  York  University 

Dept.  of  Computer  Science 

Courant  Institute  of  Mathematical  Sciences 

251  Mercer  Street 

New  York,  New  York   10012 


Work  on  this  paper  has  been  supported  by  NSF  grants  DCR-84-01898  and  DCR-84-01633. 


Polynomial  Remainder  Sequences  and  Theory  of 

Subresultants 

Chiing-jen  Ho      and      Chee  Keng  Yap  * 

Courant  Institute  of  Mathematical  Sciences 
New  York  University 

251  Mercer  Street 
New  York,  NY  10012 

August  31,  1987 


1      Introduction 

There  are  many  basic  applications  where  we  want  to  compute  the  GCD  of  a  pair  of  polynomials, 
together  with  the  eissociated  sequence  of  remainders.  These  include  Diophantine  equations,  Sturm 
sequences,  elimination  theory,  algebraic  cell  decompositions  [^'^85].  There  is  an  obvious  way  to  apply 
Euclid's  algorithm  for  integers  to  a  pair  p(x),  q{x)  of  polynomials  in  F[x]  for  any  field  F:  compute 
the  sequence 

PO,Pl>  ••  •  ,Ph,Ph  +  l  (1) 

where  po  =  p,  pi  —  q,  ph+i  =  0  and  pi^i  is  the  remainder  of  p,_i  divided  by  pi, 

p,+  i  =  rem(pi_i,pi),    i  =  l,2,...,h. 


'Thjs  work  is  supported  by  NSF  grants  DCR-84-01898  and  DCR-84-01633. 


2  1     INTRODUCTION 

Example  1 

po    =     x^  +  2x'*  +  3x2  -  X  +  2 
pi     =     3x^-1  +  2 

,       20x       16 
P.     =      3x2-— +y 

275x       166 

^^    ~     ~  243   "^  243 
115668 

Da       =         

^  75625 

P5    =     0 

However,  if  p{x),q{x)  G  I[x]  where  7  is  an  integral  domain,  the  remainder  is  in  general  not  well- 
defined.  To  handle  this,  one  defines  the  pseudo-remainder  oi p(x)  divided  by  p'(x)  to  be  the  following 
polynomial  prem(p,  p'):  if  deg(p)  <  deg(p')  orp'  =  0  then  prem(p,  p')  =  p.  Otherwise,  prein(p,p') 
is  the  polynomial  r(x)  G  I[x]  such  that  there  exists  ^{x)  €  7[x], 

r(x)  +  p'ix)q(x)  =  b'+'pix),        deg(r)  <  deg(p')  (2) 

where  6  =  lead(p'(x))  is  the  leading  coefficient  of  p'  and  d  =  deg(p(x))  —  deg(p'(x)).  Thus  by 
definition, 

prem(p, p')  =  rein(6  "'"^p,p'). 

Exercise  1  Use  long  division  to  compute  prein(p,  p')  where  p(x)  =  x®  +  2x''  +  Sx^  —  i  -|-  2  and 
p'(x)  =  3x2  -  X  +  2. 

Exercise  2  Show  that  prem(p,  p')  is  well-defined  for  all  p,p'  £  I[x]. 

Definition.  Let  p(x),5(i)  6  I[x].  Then  p  and  q  are  similar,  denoted  p  ~  g,  if  there  exists  a,b  £  I 
such  that 

ap(x)  =  bq(x) 

A  special  case  of  this  is  where  a  and  6  are  units,  in  which  case  p  and  5  aje  seiid  to  be  associates  of 
each  other. 


Definition.  Let  p{x),q(T)  be  non-zero  polynomials  in  I[x]  with  deg(p)  >  deg(7).  Then  a  sequence 
(po>  Pi  >  ■  •  ■ .  Ph)  (h  >  1)  is  called  a  polynomial  remainder  sequence  (PRS)  of  p  and  5  if  po  =  Pi  Pi  =  9 
and  for  i  =  2,3, . . .  ,h, 

p,  ~  prem(p,_2,p,-i) 

and 

prem(p/,_i,ph)  =  0. 

Following  Collins  [C66],  we  call  a  PRS  (po.Pi,  ■  ■  •  ,Ph)  the  Euclidean  PRS  if  p,  =  prem(p<_2,Pi_i) 
for  I  =  2, ... ,  h.  The  difficulty  with  the  Euclidean  PRS  is  its  unmanageable  coefficient  growth.  Thus 
the  main  issue  in  this  area  of  research  is  to  discover  efficient  methods  of  computing  polynomial 
remainder  sequences. 

Example  2   (Knuth's  example)  Displaying  only  coefficients,  the  following  shows  an  Euclidean  PRS: 

x^     x''     x^     x^        j^      x^  j^ X 1 

10       10-3-38  2  -5 

3       0  5        0-4  -9  21 

-15        0  3  0  -9 

15795  30375  -59535 

1254542875143750  -1654608338437500 

12593338795500743100931141992187500 

Remarks.   This  report  is  extracted  from  our  lecture  notes  [Y85]  (which  accounts  for  the  presence 

of  "Exercises").     It,  in  turn,   was  an  attempt  to  expound  the  theory  of  subresultants  along  the 

lines  suggested  by  Loos  [BCL83].    We  decided  to  publish  these  notes  in  view  of  the  fundamental 

importance  of  subresultants  and  the  fact  that  the  theory  is  not  easily  accessible.   Furthermore,  sis 

we  shall  point  out,  these  were  some  major  oversights  in  Loos  [BCL83]  and  the  correct  version  is 

intricate  enough  to  deserve  careful  accounting.  The  Loos-Habicht  approach  in  [BCL83]  is  basically 

different  from  the  well-known  account  of  Collins  [C67]  and  Brown  and  Traub  [BT71].  We  feel  that 

in  some  ways,  this  alternative  approach  is  more  informative  and  "digestible". 


4  2     THE  PRTMTTIVE  PRS 

2      The  Primitive  PRS 

Let  I  be  an  Euclidean  domain.   A  polynomial  p  €  I[x]  is  said  to  be  primitive  if  its  coefficients  are 
relatively  prime.  It  is  easily  seen  that  we  can  express  any  polynomial  p  in  the  form 

p{x)  =  c-q{x),        eel  -  (3) 

where  ^(x)  is  primitive.  For  instance,  over  Z[x],  p(x)  =  — 2x^  —  4x  +  12  is  not  primitive,  and  we  can 
write  p(x)  =  c  •  (x^  +  2x  —  6),  c  =  —2. 

Lemma  1    The  factorization  (S)  ts  unique  in  the  sense  thai  if  p{x)  =  c  •  q{x)  =  c'  •  q'{x)  then 

c  =  uc'  ,       u?(x)  =  q'{x) 
where  u  is  a  unit  in  7.  In  other  words,  c  and  q{x)  are  unique  up  to  associates. 

Proof.  Exercise.  Q.E.D. 

In  order  to  get  uniqueness  in  this  lemma,  it  is  conventional  to  choose  a  distinguished  element 
among  each  equivalence  class  of  associates  in  I.  Then  we  define  the  content  of  p(x)  to  be  the  element 
cont(p(x))  G  I  and  the  primitive  part  of  p(x)  to  be  the  primitive  polynomial  prim(p(x))  such  that 

p  =  cont(p)  •  prim(p) 

and  lead(prlm(p))  is  distinguished  among  its  associates.  A  beisic  result  about  primitive  polynomials 
is  the  following. 

Lemma  2  {Gauss's  Lemma)  Let  I  be  a  unique  factorization  domain.  If  p{x),q{x)  G  I[x]  are  prim- 
itive, so  is  their  product  p{x)q{x). 

Proof  Exercise.  Q.E.D. 

Exercise  3  Prove  that  if  7  is  a  unique  factorization  domain,  so  is  7[x]. 


Exercise  4  Show  that  for  p, 9  6  I[x],  cont(GCD(p, g))  and  GCD(c6nt(p),cont(g))  are  associates, 
and  prim(GCD(p, 5))  and  GCD(prim(p),prim(g))  are  associates. 

We  now  define  another  natural  family  of  PRS's.  A  PRS  (po, . . . ,  pjt)  is  a  Primitive  PRS  if  for  each 

i=  l,...,k  -I, 

prera(p,_i,p,)  ■      ,  ,  ^^ 

P.+i  =  -, -, 7T  =  prim(prem  p,_i,p,)  . 

cont(prem(p,_i,p,)) 

Exercise  5  Compute  the  primitive  PRS  for  Knuth's  example  above.  (Check:  the  final  constant 
polynomial  in  the  primitive  PRS  should  be  143193869.) 

Clearly,  this  PRS  has  the  minimal  coefficient  growth  among  all  PRS.  The  need  to  compute  the 
content  of  a  polynomial  (=  the  GCD  of  all  the  coefficients  of  the  polynomial)  could  present  difficulty. 
The  coefficient  growth  is  only  linear  in  the  univariate  case  and  this  may  be  a  reasonable  algorithm  if 
computing  GCDs  over  the  coefficient  domain  I  is  inexpensive.  However,  in  the  case  of  multivariate 
polynomials,  this  assumption  is  no  longer  true. 

3      The  Subresultant  PRS 

The  primitive  PRS  is  an  example  of  a  family  of  PRS  which  has  the  following  form:  Let 

^=(/?i,/?2,...,/?h-i),      h>l  (4) 

be  a  (possibly  empty)  sequence  of  elements  in  /.    A  PRS  {po,  ■  ■  ■  ,Pii)  is  said  to  be  based  on  the 

sequence  (4)  if  for  each  i  =  2, . . . ,  /», 

prem(p,_2,p,_i) 
Pi-i 

By  specifying  different  sequences  of  ^,  we  get  different  famiUes  of  PRS.  The  computational  properties 

of  the  corresponding  PRS  algorithm  depends  on  how  easily  the  j3i's  can  be  computed.  We  now  present 

what  is  essentially  the  best  PRS  algorithm  of  this  kind.  It  was  discovered  by  Collins  [C67]  around 

1967  with  the  present  form  due  to  Brown  [B71]. 


6  4     SUBRESULTANTS 

Definition.  Let  (po.Pi,  ■  ■  ■ ,  Ph)  be  a  PRS.  For  each  i,  let  6^  =  deg(p,')-deg(p,+i)  and  c^  =  lead(p,). 
Then  the  PRS  is  called  a  subresuliani  PRS  if  it  based  on  the  following  P  sequence: 

and  for  i  =  1 , . . . ,  /i  —  2 , 

P.^,  =  {-!)'+'•  ii'.Y'Ci  (5) 

where  the  auxiliary  sequence  (V'l ,  i/'2 .  ■  •  • .  V'ft)  is  defined  by 

^-0  =  1 

*=*-'fer  =  w;^'  •=' '  (^> 

It  is  quite  easy  to  implement  the  subresultant  PRS  algorithm.  On  the  other  hand,  it  is  not  obvious 
why  the  algorithm  works.  More  specifically,  we  have  to  show  that  each  pi  is  a  divisor  of  (the  content 
of)  prem(p,_i,  p,_2)-  Furthermore,  it  turns  out  that  the  ipi's  are  elements  of  the  coefficient  domain 
/,  in  fact  ,  they  are  determinants  of  certain  coefficients  of  po  and  pi;  this  fact  was  not  known  in 
the  early  papers  on  the  subresultant  PRS  algorithm.  These  facts  about  /?,  and  Vi  will  be  proved  in 
Section  7. 

Exercise  6  Produce  an  algorithm  that  implements  the  subresultant  PRS  algorithm.  You  may  use 
the  fact  that  the  i/'.'s  are  in  /. 

We  next  investigate  the  theory  of  subresultants  which  will  reveal  the  structure  of  PRS's  in  general 
and  of  the  subresultant  PRS's  in  particular. 

4      Subresultants 

One  of  the  first  to  realize  that  polynomial  remainder  sequences  are  intimately  connected  with  de- 
terminants was  Collins  [C66].     Our  notations  follow  Loos  [BCL83].     We  begin  by  showing  how 


pseudo-remainders  can  he  expressed  using  determinants.  For  any  sqOare  matrix  M ,  we  let  dot(A/) 
denote  its  determinant. 

Definition.    Let  A/  be  a  m  x  n  matrix  where  m  <  n.    Then  the  determinant  polynomial  of  M , 
detpol(M),  is 

det(M^)i"-'"  +  det(A/„+i)x"-'"-i  + +det{M„-i)x  +  det(Af„) 

where  Mi  (i  >  m)  is  the  square  submatrix  of  A/  consisting  of  the  first  m—  1  columns  and  ith  column 
of  A/.  We  call  det(A/m)  the  nominal  leading  coefficient  of  detpo\{M). 

Note  that  the  (actual)  degree  of  detpol(A/)  is  at  most  n  —  m  since  the  nominal  leading  coefficient 
may  vanish.  Next  we  introduce  a  convenient  notation.  Let  Ai{x)  G  I[x]  (i  =  1, . . . ,  m)  be  polynomi- 
ab  and  let  n  =  1  +  niaxi<,<mdeg(>l,).  Then  niat(^i,  A2,  . .  . ,  Am)  denotes  the  m  x  n  matrix  whose 
(:',  n  —  j)th  entry  is  the  coefficient  of  x^  in  Ai{x).  We  shall  write  detpol(>li,  A2,  ■ .  ■ ,  Am)  as  short 
hand  for  detpo\{xnBt{Ai, ...,  Am))-  The  basic  connection  between  determinants  and  polynomial 
division  is  revealed  in  the  following  lemma. 

Lemma  3   Let  A,  B  £  I[x],  deg(>l)  =  m>n-  deg(S).   Then 

detpol(x'"-"5,  x'"-"-iS,  ...,xB,B,A)  =  prem(>l,  B). 


Proof.  Write  B{x)  =  Yl"^^  bix'  and  let  A{x)  =  E^o«.^'  ^here  6  =  fc„  =  lead(B),  Let 

M  =   mat(x'"-"B,i'"-"-i5,...,xB,5,6'"-"+iyl). 

/  f>„  b„_i  ...  ...  60 

bn  ...  ...  fcl         60 

f>n  fro 

V  i'"-"+'a„     6'"-"+ia^_i      ...     b"'-"+^a„     6'"-"+iao  / 

Applying  Gaussian  elimination  that  corresponds  to  the  long-division  of  6'^~""'"M  by  B,  we  eventually 


4     SUBRESULTANTS 


reduce  M  to 

/  bn     b„_i      60 


M' 


bn  fcl  fci 


0 


bn  bo 

\  Cn-l        Oo    / 

where  the  bottom  row  corresponds  to  prein(j4,  B)  =  J2"^q  Cix'.  Here  n  —  1  is  only  a  nominal  degree 
of  prem(>l,5)  since  any  of  the  Cj 's  may  vanish.  Clearly 

detpol(M)  =  detpol(M') 

by  basic  properties  of  determinants.  We  observe  that 

detpol(A'/')  =  6"-"+^prem(^,  B) 

since  M'  is  upper  triangular.  On  the  other  hand, 

detpol(M )  -  fc'"-"+Metpol(x'"-"5,  ...,B,A). 

These  three  equations  imply  the  lemma.  Q.E.D. 

Definition.    Let  A,  B  £  I[x]  with  deg(A)  =  m  >  0  and  deg(B)  =  n  >  0.    For  any  i,  0  <  i  < 

mrn{m,n},  define  the  ith  subresultant  of  A  and  B  to  be 

sres,(vl,  B)  =  detpol(z"-'-^^,  x^-'-^A, . . . ,  A,  x'^-'-^B,  x'^-'-^B,  ...,B). 

The  ith  principal  subresuliani  coefficient  of  A  and  B,  denoted  psc^{A,  B),  is  defined  to  be  the 
nominal  leading  coefficient  of  sresi(j4,  B). 

Observe  that  inat{x"-'-'^  A,x"-'-^A, . . . ,  A,x'^-'-'^B,x'"-'-^B, . . . ,  B)  has  m  +  n  -  2i  rows 
and  m  +  n  —  i  columns  and  deg(sresi(j4,  B))  is  at  most  i.  In  particular,  sreso(^,  B)  £  /. 

For  convenience,  we  extend  the  definition  of  subresultant  sresj(^,  B)  to  the  cases  inin{m,  n}  < 
»  <  max{m,n}  —  1: 

sresi(/i,  B)  =  0       for  min{m,  n}  <  i  <  max{m,  n}  —  1 


sres,(>i,B)  =  lead(X)l"'-"l-^Y     ' 

where  /  =  nim{m,  n]  and  X  =  A  if  m  >  n;  B  i{  n  >  m. 
The  sequence 

(•S'a'  +  l.  Suj,  ■  ■  ■  ,  Si,  So) 

where  w  +  1  —  inax{m,  n],  S^^i+i  =  A,  S^,  =  B  and  5,  =  sresi(>l,  B)  for  i  =  ly  —  1, . . . ,  1,0  is  called 
the  subresuhani  chain  of  A  and  B.  It  is  convenient  to  extend  the  definition  of  principal  subresultant 
coefficient  psc,(>4,  B)  to  the  cases  I  <  i  <  w  +  I: 

psc-{A,  B)  =  lead(S,)  =  0      foi  I  <  i  <  w 

psc„,(^,  5)  =  lead(5u,)  =  lead(5)    and    psc„,^.i(A,  B)  =  1. 

Note  that  psc^^.j(yl,  B)  is  not  defined  as  lead(>l)  =  lead(S,^+i(>l,  B))  as  one  might  have  expected. 
We  remark  that  5o  is  the  classical  "resultant"  of  A  and  B  and  it  is  well-known  that  Sq  =  0 
precisely  when  A  and  B  have  a  common  non-constant  factor,  i.e.,  when  deg(GCD(y4,  S))  >  0.  The 
subresultants  generalizes  the  classical  resultant  in  that  if  Sq,  5i, . . . ,  Si_i  all  vanish  but  Si  does  not 
then  deg(GCD(>l,  B))  =  i.  These  remarks  will  be  consequences  of  theorems  to  be  shown. 

Exercise  7  1.  Let  deg(^(x))  =  m,  deg(5(x))  =  n  and  m>n>k.  Show  that  deg(GCD(>i,  B)) 
>  t  -1-  1  iff  there  exists  polynomials  Uk{x)  and  Vk{x),  not  both  zero,  deg{Uk)  <  n  —  k  —  I  and 
deg(Vlt)  <  m  —  k  —  1,  such  that 

A(x)Ut{x)-Bix)Vt{x)  =  0.  (7) 

Hint:  use  the  fact  that  A,  B  are  in  a  unique  factorization  domain. 

2.  Suppose  that  we  already  know  that  deg(GCD(>l,  B))  >  ifc.  Show  that  if  there  exists  t/t,  Vjt  as 
above,  and  Ci(x)  of  degree  at  most  ifc  —  1  such  that 

A{x)Uu{x)  -  B{x)Vt{x)  =  C{x)  (8) 


10  5     SUBRESULTANT  CHAINS  WITH  INDETERMINATE  COEFFICIENTS 

then  (7)  has  a  non-trivial  solution.  (The  reverse  implication  ii  trivial.) 

3.  Show  that  (8)   is  equivalent  to  psci^{A,B)  -  0.     Conclude  that  deg(GCD(yl,  5))   =  t  iff 
pscj,(^,  B)^  0  and  for  all  !  =  0, . . . ,  jt  -  1,  psc,(yl,  B)  =  0. 

5      Subresultant  chains  with  indeterminate  coefficients 

We  next  follow  Loos  [BCL83]  in  developing  the  basic  results  of  subresultants  in  terms  of  polynomials 
with  indeterminate  coefficients.  In  other  words,  for  this  section  we  assume  that  the  input  polynomials 

m  n 

^=^a.x-,        5  =  ^6.:r'  (9) 

t=0  1=0 

come  from  the  ring 

Z[z;  am,  ■  ■  .,ao,b„,.. . ,  6o] 

where  a,  and  6,  are  indeterminates.  After  obtaining  the  properties  of  subresultants  in  this  setting,  we 
can  interpret  the  indeterminate  coefficients  of  yl  and  B  as  values  in  I,  thereby  observing  the  behavior 
of  the  subresultant  chain  in  I[x].  This  approach  has  several  advantages  as  noted  by  Loos  [BCL83]. 
In  particular  it  separates  the  sparsity  effects  (when  consecutive  chain  members  have  degrees  differing 
by  more  than  one)  from  the  similarity  relations  among  various  polynomials  defined  from  the  chain. 
In  this  section,  we  assume  that  the  polynomials  A  and  B  in  (9)  satisfies  m  =  n  +  1.  We  will  see 
that  this  is  without  loss  of  generality  in  our  applications. 

Lemma  4  Assume  that  the  polynomials  A  and  B  in  (9)  satisfies  m  =  n  +  I.   Then 

(i)  sreSn--i{A,  B)  =  pTein{A,B). 

(ii)   Fork  ^0,...,n-2,  bl^"~''~'^\reSk{A,  B)  =  sresi(5,prem(^,  B)). 

Proof.  Part  (i)  is  just  a  restatement  of  the  previous  lemma  since 

sres„_i(^,5)   =   detpo\{A,xB,  B) 

=  detpol{x  B,B,  A)  =  prein{  A,  B). 


11 


To  see  part  (ii), 

sresi(^,  S)     =     detpol(x"-^-''A,...,A,T"~''B,...,B) 

=     detpol(x"-'-HlA, ...,  blA,  x"-*B B)  ■  i"^'"-*) 

=     dctpol(i"-'-*prem(yi,  B),...,  prem(>l,  5),  i"-*B, . . . ,  B)  ■  6;:^^""*^ 

(by  subtracting  appropriate  multiples  of  S-rows  from  yi-rows.) 
=     detpol(z"-*5, . . . ,  5,  x"-i-'=prem(A,  S), . . . ,  prem(>l,  B))  ■  6;^^""*^ 

(by  transposing  rows) 
=     detpol(i"-*-2S, .  . . ,  B,  x"-i-*^prem(A,  B), . . . ,  prem(^,  B))  ■  6;^'""*^"^^ 

(using  the  fact   that   ib  <  rz  —  2  and  deg(prein(j4,B))  =  n  —  1) 
=     sresi(B,prem(>l,B))  ■  bn  ~  '. 

Q.E.D. 

We  note  that  Lemma  4  (i)  justifies  viewing  subresultants  as  a  generalization  of  pseudo-remainder. 

Theorem  5   (  Habichi's  theorem  )  Let  Sn+i-,  Sn, . . . ,  Sq  be  the  subresultant  chain  of  A  and  B.   Let 
Rj  (j  =  0, . . .  ,n  +  1)  be  the  jth  principal  subresvHant  coefficient  of  A,B .   Then  for  all  j  =  1, . . . ,  n, 

Rf^;^^St  =  sresi(S;  +  i ,S,)  ,     k  =  0, .. .  J  -  I  '     (10) 

R'j  +  iSj-i  =  prem(5j  +  i,5j).  (11) 

Proof     We  first  observe  that  (11)  follows  from  (10)  since  substituting  j  —  I  for  k  in  (10)  gives 
R]+iSj-i  =  sreSj-i{Sj+i,Sj)  =  prem(5j+i,  Sj)  by  the  last  lemma.    It  remains  to  show  equation 
(10).  Basis:  If  j  =  n  then  equation  (10)  follows  from  the  definitions.  Induction:  let  j  <  n,  assuming 
the  lemma  is  true  for  j  +  1.  Let  k  <  j  —  1.  By  the  induction  hypothesis 
^20-i+i)^^     _     sresi(Sj+2,5'j+i) 

=     -RJ+i        sreSi(Sj.|.i,  prem(Sj+2,5j+i))       (by  previous   lemma  (ii)) 

2. 


^'-'\rest{S,+:,R]^,Sj)      (by(ll)) 
(^■-'^detpoI(xJ-*-iS;+i, . . . ,  Sj+,,R]^^x^-'S, ,...,  i?j%25j) 

This  completes  our  proof  of  (10).  Q.E.D. 
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6      Specialization  of  subresultant  chains 

We  now  want  to  substitute  the  indeterminate  coeflRcients  a,,  bj  in  Z[r;  On+i, . . . ,  oq,  fc„,  •  •  • ,  io]  with 
specific  values  in  the  domain  7.  This  is  called  specializing  the  indeterminate  coefficients  and  it 
induces  a  natural  homomorphism  from  Z[x;  a„+i, . . .  ,ao,bn,  ■  ■  ■  ,bo]  to  I[x].  Let  (t  denote  any 
such  homomorphism.  Let  S^  (resp.  Rt)  denote  the  kth  subresultant  (resp.  principal  subresultant 
coefficient)  of  A  and  B  regarded  as  elements  of  Z[x;  On+i, . . .  ,ao,b„, . . . ,  bo].  Let  A,  B,  St  and  Rt 
denote  the  image  of  A,  B,  Sk  and  Rt  under  a.  In  general,  we  write  a  over  an  expression  to  denote 


the  image  of  that  expression  under  <t,  e.g.  sres,(Sj+i,  Sj).  It  is  important  to  realize  that  this  is  not 
the  same  as  sresi(Sj+i,  5j).  This  fact  seems  to  have  been  overlooked  by  Loos  [BCL83]  therefore 
rendering  his  proofs  invalid.  Our  goal  is  to  investigate  the  specialized  chain 

(5n+l  ,  ■  •  •  ,Sq). 

Let  dk  and  k  denote  the  actual  and  nominal  degree  of  St  respectively,  with  the  usual  convention  that 
deg(O)  =  —GO.  A  subresultant  Sjt  is  said  to  be  defective  of  degree  r  \i  dt  =  r  <  k;  otherwise  St  is 
said  to  be  regular.  A  chain  is  said  to  be  defective  if  any  of  its  members  is  defective;  otherwise  it  is 
regular.  Thus  for  a  defective  chain,  the  sequence  ((f„+i,d„, . . .  ,do)  contains  gaps  in  the  sense  that 
some  integer  in  the  set  {0, .  . .  ,  n}  is  missing.  The  growth  of  the  content  of  polynomials  in  any  PRS 
is  closely  linked  to  the  presence  of  gaps.  Note  that  lead(Si)  /  R^  if  St  is  defective  or  Jk  =  n  +  1. 
Definition.  Let  A,B  £  I[x]  with  deg{A)  =  m  >  deg(S)  =  n  >  0.  For  any  j,  0  <  t  <  m  —  1,  define 
the  ith  ■pseudo  subresultant  of  A  and  B  to  be 

psres,(^,  B)  =  detpo\{x"'-'-'^ A, x'"— ^yi, ....  A,x"'-'-'^B,  i"'— ^5, ...  ,5). 

Note  that  if  Sjj^i  is  regular,  then  for  0  <  it  <  j 

psreSi(5j  +  i,S,)  =  detpol(x>-i-*Sj+i, . . .  ,S;  +  i,x^-*S,, . . .  ,5,)  =  sresi(Sj+i,Sj). 
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2  j  +  1  -  ib 


r  —  t  +  1  rows  of  Sj+\ 
(j  -  r  -  1)  0  V 


j  —  k  rows  of  Sj+i 


j  —  k  +  I  rows  of  Sj 


r  +  2 


Figure  1:  The  matrix  associated  to  psresi^{Sj+i,Sj) 

If  both  Sn+\  and  Sn  are  defective,  then  it  is  easily  seen  that  the  rest  of  the  chain  is  0.  Henceforth, 
we  assume  that  at  least  one  of  Sn+i  or  5„  will  be  regular.  Without  lost  of  generality,  we  may  aissume 
that  Sn+i  is  regular.  The  following  lemma  is  the  specialization  of  Lemma  4: 

Lemma  6  Lei  Sn+i,Sn, . . . ,  So  be  the  pseudo  subresultani  chain  of  A  and  B.  Lei  (j  be  the  actual 
leading  coefficieni  of  Sj.  For  j  =  l,...,n,  ifSj+i  is  regular  and  Sj  has  degree  r+l  {—oo  <r+l  <  j) 
then 


(i)  psres,(S,  +  i,S,)  =^j;[-^(-lV-'+iprem(S;  +  i,S;). 


c.\  _  <rJ-'--lf(i-'--l)e, 


(ii)  psres,+  i(5,  +  i,S,)  =  ^^^r'^''^     'S, 


(iii)   For  k  =  0, . . .  ,r  —  1, 


psreSi(Sj+i,Sj) 
=     psres,(5,,prem(S,+i,S,)).^j;r'-(-l)^^+'-'^^^-*+'^-^7^'"'"''^^'''^ 
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Proof,  (i)  By  Figure  1,  we  have 

psres^(Sj  +  i,Sj)   =  fj~[~Metpol(5^  +  i,x^-''5j, . . .  ,5_,) 

=  .ej;r'(-iy-'+'prem(Sj  +  i,Sj)     (by  Lemma  3) 
Also  from  Figure  1,  we  can  ecisily  get  (ii).  To  see  part  (iii), 
psreSi(Sj+i,5j) 
=     detpol(z-'-^-*5j  +  i,  ...,Sj  +  i,x'-''Sj,...,Sj) 
=     detpol(z>-i-*S,  +  i, . . . ,  x^-'^+^Sj+i,  x^-*^j-^+^5,+i, . . . , ^j-^+^5,+1,  x^-*S,-, . . . , 5,) 

(j-r+l){r+l-k) 

=     detpol(x-''~^~''5j  +  i, ..  .,z'"~'=  +  ^5j+i,i'"~*prem(Sj+i,5_,), .. .  ,prem{Sj+i,Sj),x'-^Sj, . . .  ,Sj) 

c-U-r+i)(r+i-k) 

=     detpol(r-'-i-*Sj  +  i x'-*+iS^  +  i,x^'-*S_,-,...,i'-*S;,x'-*-^Sj,...,Sj,x'-'=prem(Sj+i,5j),. 

prem(5j  +  i,5,))-(-l)(^+i-^)0-'^+i).^-^^-^+''''+'-*) 
=     detpol(x'"~*'~^Sj, . . .  ,Sj.  x''~*prein(Sj4.i,5j), . . .  ,prein(5j  +  i,  Sj)) 

_^j-r-l       fj-r+l      ,_^yr+l-i)(;-i+l)  .  £-(; -'"+1  )(r+l-l) 

=     psres,(5,,prem(.S  +  i,5,)).^j;r^-^i-'-+^-(-l)(^+^-')(^-'+^)-^-^^-^^^^^^^^-'^ 
Q.E.D. 

Theorem  7  (  Specialized  Habichi  's  theorem  )  Lei  Sn+i,  S„, . . . ,  So  be  the  -pseudo  suhresuliant  chain 
of  A  and  B.  Lei  Rj  (j  =  0,...,n  — 1)  be  ihe  nominal  leading  coefficient  of  psreSj{A,  B).  Let 
R„  =  lead(5)  and  Rn+i  —  1.   i^or  j  =  1, .  . .  ,n,  tf  Sj+i  is  regular  and  deg(Sj)  <  j  then 

-R^r^'-^i  =  psreSi(S,+i,  5,)  ,     k=0,...,j-l  (12) 

Proof.    We  prove  the  theorem  by  induction.    Basis:    If  j  =  n  then  equation  (12)  follows  from  the 
definitions.  Induction:  assuming  the  lemma  is  true  for  j+l  and  Sj  has  degree  r+l  (— oo  <  r+1  <  j). 
We  want  to  show  that  .Rr+T    -^it  =  psrest(Sr+i,Sr)  for  ib  <  r  —  1.    By  the  induction  hypothesis 
Rj.^y      Sic  =  psresjt(S'j4.i ,  5j).  We  have  two  cases  to  consider: 
Case  (i)  :  j  —  n,  then  jR;+i  =  1.  Thus,  by  the  previous  lemma  (ii),  we  have 

Sr+i  =  psres,+  ,(S,+i,5;)  =  ^^'^'^'^'''''^3^. 
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Hence, 

^;  =  U^;  +  i^;)"'^"'"'^-^r+i     and     ^.+  i  =  ^j^r^^j"^  (13) 

Also  by  the  previous  lemma  (i),  we  have 

Sr  =  psres,(S,+i,S,)  =  ^j;[-^(-l)'-^+iprem(S;+i,S,). 


Thus, 


(-^;  +  i)-'^-'-'^5.  =  prem(S,  +  ,,S,)  (14) 


Hence, 

Sk 
=     psresj.(5'j^.i ,  5j)       (by  inductive  hypothesis) 

=     psres,(S„prem(5,>i,5,))  .^j;r^  •  (-l)(^+^-*)(^-'=+^)  •  ^;^^-^+^^(^-^) 
=     detpol(x'"~*~^Sj, . . . ,  5j,  x'""*prein(Sj4.i,  Sj), . . . ,  prein(5j+i,  5j)) 

.^j-[-i.(_l)(r+i-i)(;-J:+i).^-(j-+i)(--'^)       (by  previous  lemma  (iii)) 
=     detpol(x'-^-i(^,+i^;)-(^-^-i)S,+  i,...,(^,+i^,)-0-^-i)5.+  i,x^-*(-^,  +  i)-(>-^-i)5„..., 

(-^,+i)-(^-'-^)S.).^j;[-^.(-l)(^+i-*)(>-^+i).^-^^'-^+^)(^-*^    (by  (13)  and  (14)) 
=     detpol(x^-*-^S.+  i Sr+ux'-'Sr, . . .  ,S.)  •  (^;+i  ■  ^,.)-0  — i)(-i)  .  (-^^.^i)-(J  — i)(-i+i) 

.^^i-'"-! .  (-nCr+i-tja-t+i) .  ^-o-'-+i)('--'=) 

-  r^=T-^c    l'9  ^^      c-2(;-r-l)(r-J:)         _(;-r-l)(r-i) 
.(_l)(j-'-+l)('--i  +  l)   .  (_l)(r+l-i)(;-k+l)   .  ^-U-r+l^r-k) 

=     psreSi(S.+  ,,5,)-^-;p-^-''(^-').(-l)(^+i-*X2^-2-+2+-'^).e-'^^-'^^'"-'^ 

-  T.«,.pe    (^  ^    \      ,-2U-r-l){r-k)         -2(;-r)(r-i) 

=     psTest{Sr+,,Sr)-R;+['~'^     (by  (13)) 
Case  (ii)  :  j  ^  n,  then  ^j+i  =  ^>+i-  Thus,  by  the  previous  lemma  (ii),  we  have 

Hence, 

5,  =  ^;:;r^^-(^-^-^)s.+  i    and    ^+1  -  ^-|r^-^)^j-^  (15) 
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Also  by  the  previous  lemma  (i),  we  have 

Rf^r^S,  -  psres,(5,+i,S,)  =  R^-[-\-iy-'+'pr^Tn{S,+uS,). 


Thus, 


Hence, 


i-Rj  +  ^y-'+'Sr  =  prem(5j+i,S,)  (16) 


=     psres^(5y+i,  5j)       (by  inductive  hypothesis) 

=     psres,(5,,prem(S,>i,5,)).^j;r^-(-l)('-+^-^)(^-^+^)-^7"-^+^^^^-'^ 
=     detpol{x''~''~^Sj, ..  .,S_,-,2:''~*prem(5j+i,5j),.. . ,  prem(Sj+i,  5^)) 

■R^+r'  ■  (-l)(^+i-»^K^-^+»)  •^7^^-'+''('-')       (by  previous   lemma  (iii)) 
=     detpol(x^-»^-i^j;r '^-(^-^-^^S.+  i, . . . ,  R]l[-'i-^^-'-'^Sr+i,x-'{-Rj+ry-^+'Sr, .... 

(-^j+i)>-^+i5,)    .Rj;p'-(-l)(^+i-*)(^-'=+i)-,e7''~'^'^^'"'^    (by  (15)  and  (16)) 
=     detpol(x-*-iS.+i, . . . ,  S.+  i,  r-^S., . . . , S.)  •  R^^-r'^^^'^  ■  ^-(^-^-iX^-*)  .  {-R,^,)U-r^iXr-t+i) 

.Ri"'-'^ .  (_i)('-+i-'fc)(j-i+i) .  ^-(j-'-+i)('-*=) 

=     psres,(S.+  :,S.)  •  ^a-;-i)(-^+i)  .  ^-(i-r-iKr-.)  .  ^o--r+i)(.-.+i) 

.(_l)(;-'-+n(r-i+I)  .  (_i)(r+l-J:)(j-i+l)  .  ^-(i-'-+l)(r- 1) 

=     psres,(S,+  i,S,)-^-VT'^^'"'^'^-(-l)''^'"'^'''"''^'^'"*^-'^7'^'"'^^''"'^ 
=     psres,(S.+  „5.)  •  fi^O-^-DC.-.)  .  ^.O-^)  .  ^-2(;-r)(r-.) 

=     psres,(S.+  i,5,).^;;(;-*).ftj(.V    (by  (15)) 
This  completes  our  proof  of  (12).  Q.E.D. 

Theorem  8  {Pseudo  Subresultant  Theorem)  Let  (Sn+i,  ■  ■  ■  ,§0)  be  a  specialized  chain  in  I[x].   For 
j  =  1, . . . ,  n,  t/Sj+i  15  regular  and  Sj  has  degree  r  (—00  <  r  <  j)  then 

1.  j  =  n, 

5j_i=5,_2  =  ---  =  S,+  i=0  (17) 

5r  =  lead(5j)-'-'"lead(Sj+i )■'■-'" 5j,      0  <  r  <  j  (18) 
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5,_i  =(-lead(5,  +  i)V-^prem(Sj  +  i,5,),  •    l<r<;  '      (19) 


2.  j  <  n, 


Sj.,  =  Sj.2  =  ---  =  Sr+i^0  (20) 

{Rj  +  ,y-'Sr  =  \eBLd{S,y-'Sj,      0<r<j  (21) 

{-Rj+,y-'-+^Sr-,=prexniSj+uSj),      l<r<j  (22) 
Proof.  From  (12)  we  see  that  for  k  =  r  —  1,  r, . . . ,  j  —  1, 

>  RfS'^S,=psTes,{S,  +  ,,S,).  (23) 

The  right-hand  side  of  (23)  is  equal  to  detpol(Af )  where  M  has  the  form  as  shown  in  Figure  1. 
(Note:  In  Figure  1,  the  degree  of  Sj  is  r+  1.)  Observe  that  S,  has  j  —  r  leading  terms  vanishing,  and 
there  are  j  —  k  rows  of  S^+i  in  the  matrix  M.  If  j  —  r  >  j  —  k,  i.e.  k  >  r,  then  it  is  clear  from  the 
picture  of  Af  that  M  will  be  upper  diagonal  (i.e.  has  only  zeroes  below  its  main  diagonal).  Observe 
that  the  coefTicients  in  detpol(A/)  are  all  determinants  of  upper  triangular  square  matrices.  We 
consider  three  caises. 

Ca^e  k  =  r  and  j  =  n:  Then  .R;  +  i  =  1  and  detpol(M)  evaluates  to  lead(Sj+i)-'~'"lead(Sj)-'~''Sj. 
This  shows  (18). 

Case  k  =  r  a.nd  j  <  n:  Then  detpol(M)  evaluates  to  .Rj~[Iead(S_, )■'-'■  5,.  This  shows  (21). 

Case  k  >  r:  Each  diagonal  of  the  square  matrices  determining  coefficients  of  detpol(M)  has  at  least 
one  zero  element,  hence  each  coefficient  is  zero.  This  shows  (17)  and  (20). 

Now  suppose  j  —  r  =  j— k—  i,  i.e.    k  —  r  —  1.    Then,  on  expanding  the  first  j  —  r  rows  of 

detpol(M),  we  obtain 

detpol(M)   =  lead(Sj+i)^-'"detpol(Sj+i,i-'-'"+iSj,2^-'"5_,, . .  .,Sj) 

=  (-lead(S;+i))>-^detpol(i>-^+iS;,...,S,,S;+i)  (24) 

=  (-lead(5^  +  i))-'-'"prem(Sj  +  i,Sj). 
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If  j  =   n,   then  Rj  +  i    =   1.     This,  combined  with  (23)   and  (24),   ImpHes  (19).     If  j    <   n,   then 
^_,+  i  =  lead(5j  +  i).  This,  combined  with  (23)  and  (24),  implies  (22).  Q.E.D. 

Let  (Sn+i,S„, . . .  ,So)  be  the  pseudo  subresultant  chain  of  A  and  B.  Let  (S„+i, . .  .,5o)  be  the 
subresultant  chain  of  A  and  B.  Let  Rj  {j  —  0, . . .  ,n  +  1)  be  the  jth  principal  subresultant  coefficient 
of  A,  B.  Suppose  the  degree  of  S  is  /  <  n  +  1.  Then  clearly 

Sj  =  lead(S„+i)"-'5j      and      Rj  =  lead(S„+i)"-'/2j      for     0  <  j  <  n.  (25) 

Before  showing  the  subresultant  structure  theorem,  we  note  the  following  useful  consequence  of 
lenuna  3: 

Corollary  9   Lei  A,  B  £  I[x]  and  a,  0  E  I ■   Then 

pTem{aA,PB)  =  Q/?-''"''+^prem(.4,  B) 

where  j  +  I  =  deg(A)  >  r  =  deg(B). 

Theorem  10  {Subresuliani  Structure  Theorem)  Lei  {Sn+i,  ■  ■  ■  ,So)  be  the  subrtsuliant  chain  of  A 
and  B.  For  j  =  1, . . .  ,n,  tf  Sj+i  is  regular  and  Sj  has  degree  r  (— oo  <  r  <  _;')  then 

5;_,  =5j_2  =  •••  =  5,+  i  =0  (26) 

(i2_,  +  i)-'-'"5p  =  lead(S>)^-'"5;,      r>0  (27) 

i-Rj+iy-'+^'Sr-i  =  prem(5;+i,S;),      r  >  1  (28) 


Proof.  From  (20)  and  (25),  we  immediately  get  (26).  To  see  (27),  we  consider  the  following  two 
cases.  Ccise  one:  if  j  =  n,  by  definition  Rj+i  =  1  and  5j+i  =  Sj+i  and  Sj  =  Sj.  Thus,  the  right 
hand  side  of  (18)  is  equal  to  lead ( 5_,+i)-'-'" lead (5j)-'-'"5j  and  by  (25)  the  left  hand  side  of  (18)  is 


19 
equal  to  lead(5j+i)-'-'"Sr  =  lead(Sj  +  iy-'"i?j;[S,.  Hence, 

R]~[Sr  =  \ead{Sjy-'Sj       for  j  =  n,  r  >  0  (29) 

Case  two:  if  j  <  n,  then  by  (25),  (21)  can  be  evaluated  to 

(lead(S„+i)"-'/Zj+i)^-^lead(S„+i)"-'S,  =  lead(lead(S„+i)"-'5j)^-^lead(5„+i)""'5j 

(Note  that  this  equation  is  true  even  when  j  +  I  =  n  because  in  this  case  S^+i  is  regular  and  then 
n  —  I  and  thus  lead(Sn+i)"~'  =  1  so  that  we  can  add  it  anywhere  in  this  equation.)  Hence,  we 
have  when  j  <  n 

{Rj  +  ^y^Sr-  =  lead(5j)-'-'"Sj,       for  ;  <  n,r  >  0 

This  with  (29)  shows  (27). 

To  see  (28),  we  also  consider  the  following  two  cases. 
Case  one:  if  j  =  n,  by  definition  Rj+i  =  1  and  Sj+i  =  ■S'j+i  and  Sj  =  Sj.  Thus,  the  right  hand  side 
of  (19)  is  equal  to  (  — lead(Sj+i))-'~'"prem(5'^  +  i,  5^)  and  by  (25)  the  left  hand  side  of  (19)  is  equal 

to  lead(S_,  +  i)^-'-Sr_i  =  lead(Sj+iy-^ R^'l'^^Sr-i.  Hence, 

{-Rj+ry-'+-Sr-i  =  prem(5,  +  i,5j),       for  j  =  n,r  >  1  (30) 

Case  two:  if  j  <  n,  then  by  (25)  and  Corollary  9,  (22)  can  be  evaluated  to 

(-lead(5„+i)"-'iEj+i)'"'+'lead(5„+i)"-'5,_i 
=  prem(lead(S„+i)"-'5j+i,lead(S„+i)"-'Sj) 

=  lead(5„+i)""'(lead(5„+i)"-')^-^+'prem(S,+i,S,) 

(Note  that  the  above  equation  is  true  even  when  j  +  I  =  n.)  Hence,  we  have  when  j  <  n 

i-Rj+iy-^-^^Sr-i  =  prem(S,>i,5j),       for  ;  <  n,r  >  1 
This  with  (30)  shows  (28).  Q.E.D.  ' 
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Note:  This  theorem  is  essentially  the  "Subresultant  Theorem"  rn  Loos  [BCL83]  (page  122).  The 
subresultant  structure  theorem  describes  the  gap  structure  in  subresultant  chains.  In  particular,  it 
says  we  can  uniquely  divide  the  chain  into  blocks  of  (consecutive)  subresultants  such  that  if 

(5.,S,_i,...,S;+i,S;)  ,     n+1  >»■>  j>0, 

is  a  block  then 

1.  Either  j  =  0  and  all  the  subresultants  in  the  block  are  zero,  or  else  both  Si  and  Sj  are  non-zero, 

Si  ~  Sj  and  S,_i  =  5i_2  =  •  •  ■  =  5^+1  =  0.  Note  that  there  is  at  most  one  block  where  all  the 
subresultants  are  zero;  call  this  the  zero  block. 

2.  Assume  that  the  block  under  consideration  is  a  non-zero  block.  Then  Sj  is  regular.  Hence  if 
i  >  j  then  Si  is  defective.  Thus  every  defective  subresultant  Si  of  degree  r  corresponds  to  a 
unique  regular  subresultant  Sj,  r  —  j,  that  is  similar  to  it.  This  also  implies  that  there  cannot 
be  two  consecutive  non-zero  subresultants  that  are  both  defective. 

7     Relation  between  Subresultant  Chain  and  Subresultant 
PRS 

Now  we  relate  the  subresultant  PUS  p  which  is  described  in  Section  3  (equations  (5)  and  (6))  to  the 
subresultant  chain 

S  =  (S„+i,5n,  .  ■  ■ ,  So) 

where  5n+i  =  po  and  Sn  =  pi-  (Note  that  the  convention  for  subscripting  for  PRS's  in  increasing 
order  is  opposite  to  that  for  subresultant  chains.)  By  the  definition  of  PRS,  deg(po)  >  deg(pi). 
In  Section  6,  we  only  discuss  the  Ccise  deg(po)  >  deg(pi).  Now,  we  need  to  consider  the  case 
deg(po)  =  deg(pi). 

Suppose  A,  B  €:  I[x]  and  deg(yl)  =  deg(B)  =  n.    Let  (Sn+i,  5„, .. .  ,So)  be  the  subresultant 
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chain  of  >1  and  B.    Let  {Rn+i,  Rn,  ■  ■  ■ ,  Ro)  he  the  principal  subresCiltant  coefficients  of  A  and  B. 
Then 

Lemma  11  S„_i  =  — prem(j4,S). 

Proof. 

Sn-i     -     sreSn-iiA,  B) 

=     detpoliA,B)      =     -detpol{B,A)     =     -prem(^,5) 

Q.E.D. 

Let  deg(Sn_i)  =  /  and  lead(Sn)  =  ^. 
Lemma  12   For  j  =  0, ...,/,  S,  =  ^'"'sres, (Sn,  Sn-i). 
Proof. 

Si       =       STeSi{Sn+l,Sn) 

=  detpol(x"  — i5„+i, . . . ,  S„+i,  i"— iS„, . . . ,  S„) 

=  detpol(r"  — ^^S„+i , . . . ,  ^S„+i,  z"— iS„, . . . ,  S„)  •  C"" 

=  detpol(x"-'-^prem(5n+i,S„),.. . ,  prem(S„+:,  S„),  i""'~^S„, . .  .,5„)  ■  ^'-" 

=  detpol(x"— i5„_i , . .  . ,  S„_i,  x"-'-i5„, . . . ,  S„)  ■  {-0''" 

=  detpol(x"— i5„,...,S„,x"-'-i5„_i,...,S„_i)-(-0-"-(-l)<"-'^^"-'' 

=  detpol(x"  — iS„, ....  x'-'5„,  x'-'-i5„, . . . ,  S„,  x"— iS„_i, . . . ,  S„_i)  •  C'" 

=     sres.(5„,5„_i)-r"'-r-" 
=     f'-'sres,(5n,S„_i) 
Q.E.D. 

Let  {Pn,Pn-i,  ■  ■  ■ ,  Po)  be  the  subresultant  chain  of  Sn  and  Sn-i-  Let  (Cn,  C„_i, . . .  ,Co)  be  the 

principal  subresultant  coefficients  of  S„  and  Sn-i  •  Then,  from  Lemma  11,  for  i  =  0, . . . ,  /, 

Si  =  C~'Pi  and  Ri  =  e~'Ci.  (31) 

Theorem  13  Suppose  A,  B  £  I[x]  and  deg{A)  =  deg(5)  =  n.  Lei  {Sn+i , . .  .,Sq)  be  the  subresul- 
tant chain  of  A  and  B.  Suppose  deg(5n_i)  =  /.   Then, 

5;  =  lead(S„_ir-'-'5„_i  (32) 
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(-l)"-'+i^5,_i  =prem(S„,5„_i)-  (33) 

For  j  =  1, .  .  .  ,  n  —  2,  if  Sj  +  i  is  regular  and  Sj  has  degree  r  (— oo  <  r  <  ;)  then 

5,_i=Sj_2  =  ---  =  S,+  i=0  (34) 

(/?,  +  i)^-^S.  =  lead(S;)^-^S;,       r>0  (35) 

(-i?,  +  i)^-^+^5,_i  =  prem(S,+i,S,),      r  >  1  (36) 

Proof.  Because  deg(5n)  >  deg(Sn_i),  we  can  apply  Theorem  10  so  that  for  j=l,...,n—  1,  if 
Pj+i  is  regular  and  P,  has  degree  r  (— oo  <  r  <  j)  then 

P,_l=P;_2=-=Pr+l=0 

{C,  +  ,y-'Pr  =  le^d{Pjy-'Pj,      r  >  0 

(-Cj  +  i)^-^+'P,_i  =  prem(P,+i,P,),      r  >  1 

Thus,  Pi  =  lead(P„_i)"-i-'Pn_i.  Hence,  by  (31),  we  have  5,  =  lead(S„_i)"-^-'S„_i  which 
proves  (32).  Since  (-l)"-'+ip,_i  =  prem(P„,  P„_i)  ,  we  have,  by  (31),  (-l)"-'+^^5(_i  = 
prem(5n,  5n_i)  which  proves  (33). 

Also  by  (31),  we  have  for  j  =  1, ...,/—  1,  if  Sj^i  is  regular  and  Sj  has  degree  r  (— oo  <  r  <  j) 
then 

^'-(^-^)S,_i  =  ^'-(^•-2)S,_2  =  •  ••  =  ^'-('+^)5.+  :  =  0 

U'-^^  +  '^R,  +  ,y-'^'-'Sr  =  lead(^'->S,y-^^'-^5,,      r  >  0 

(-^'-(^■+^)p,  +  i)^-^+V-(^-'^5._:  =prem(^'-0>i)5,>i,e'-^5,),      r>  1 

But  ^  in  these  equations  can  be  canceled,  thus  we  get  (34)  (35)  and  (36).  Q.E.D. 

Recall  the  block  structure  of  subresultant  chains.  We  make  the  following  observations: 
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Lemtna  14        L    The  number  of  non-zero  blocks  in  the  chain  S  is  tqual  to  the  length  h  +  \  of  the 
subresullani  PES  p.    Hence  wr  can  write  S,  omitting  its  zero  bhck  if  any,  in  the  form 

S  =  So;  Si;  . . . ;  S^ 
where  St  is  the  kih  block. 
2.   For  each  k,  let  5rf(jt)  denotes  the  regular  subresultant,  clearly  0  <  d{k)  <  n  +  1,  in  the  kth  block 
Sk  is  similar  to  pi^  m  the  subresultant  PRS. 

Proof.  We  show  by  induction  on  k  that  the  zth  polynomial  in  p  is  similar  to  Sa(^i)  for  t  =  0, . . . ,  h. 

The  result  is  true  {or  k  —  I  since  h  >  2  implies  po  and  pi  are  in  separate  blocks.  If  t  >  2, 
Pk     =     prem(pt_2,pi:-i) 

~     prem(S<f(jt_2).  ^^(i.i))  (by  induction   hypothesis) 

~     ■S'(i(jt)  (from  the  subresultant  structure  theorem) 

Q.E.D. 

We  now  want  to  find  the  coefEcients  of  similarity  between  5^(4)  and  pt,  k  =  0,  ...,/i.  By 
definition,  the  first  subresultant  in  the  block  Si  is  Sd(k-i)-i  (where  we  take  d{—l)  —  n  +  2).  It  is 
convenient  to  let  e(ib)  =  (i(Jt  —  1)  —  1.  Call  (5j(o), . . . ,  Sc{h))  the  subregular  subresuliants.  Note  that 
all  the  nonzero  defective  subresultants  are  included  among  the  subregular  subresultants. 

Theorem  15   (  Correctness  of  the  subresultant  PRS  )  The  subregular  subresultants  in  S  is  related 
to  the  subresultant  PRS  p  by 

(i)   i?d(i)  =  V'i.  i  —  0,1, . .  .  ,h.  {except  that  ^&i  =  1  when  deg(po)  =  deg(pi).) 

(ii)  5e(,)  -  Pi,  i  =  0,l,...,h. 

Proof.  We  proof  both  parts  by  induction  on  i. 

Basis:  If  i  =  0  then  Rd(o)  —  Rn+i  =  1  =  V'o  by  the  definitions.  Also,  by  definition,  •5e(i)  =  S„  =  pi 
and  Sj(o)  =  5<i(o)  —  Sn+i  -  po-  Thus 

_  prem(po,Pi)  _  prem(5<i(o),  5,(1))  _  (-l)*''+^5s(2)  _ 
^-  Jr A  -       (-1)^0+1       -^«(^) 


21  7    RELATION  BETWEEN  SUDRESULTANT  CHAIN  AND  SUDRESULTANT  PRS 

If  deg(po)  =  deg(pi),  then  Sq  =  0,  so  that  ^i  =  1  and  ip2  =  c*'-  From  (32),  we  know  that 

Rj(2)  =  lead(5,)  =  cp  =  V'2. 

And  from  (33),  we  have 

(-l)"-'+^e5.(3)  =  prem(S,(i),  5e(2)) 

Hence, 

prem(S.(i),Se(2)) 

Inductive  steps:  for  J  >  1,  assuming  that  part  (i)  is  true  for  i  —  1  and  part  (ii)  is  true  for  i  and 
2  +  1,  we  want  to  show  that  (i)  hold  for  i  and  (ii)  hold  for  i  +  2.  We  rewrite  (27)  by  using  the 
notation  d{i)  and  e{i)  as  follows: 

(^d(.-i))*--'5d(.)  =  lead(S,(.))*--'5.(i)  (37) 

Thus  (V',-i)'-'-'5d(,,  =  Iead(p.)'-'-V.  and  (V'.-i)'-'-'-Rd(O  =  c'-'.  Hence, 

which  proves  part  (i).  To  see  part  (ii),  we  rewrite  (28)  as  follows: 

{-Rd(i))  '"'"^•S'c(i+2)  =  prem(Sd(,),Se(,+i))  (38) 

Then 

A+iPi+2  =  prem(p,,pi+i) 

=  prein(Se(,),Se(,^i))      (by  inductive  hypothesis.) 


"P'^^'°^ead(5.,,))^-.-i-^^^->^'<'^^)^      ('^^(''^^ 

=  — ^ — 3^prein(Sd(,),Se(,.).i))      (by  inductive  hypothesis  and  CoroUary  9.) 


=  ^^Hrr(-V'.y'^'-?<(i+2)      (by  (38)  and  part  (i)) 
Hence,  5e(,+2)  =  Pi+2,  which  proves  part  (ii).  Q.E.D 
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Corollary  16    The  il-\  's  defined  in  (6)  are  elements  of  the  integral  domain  I. 
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